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The properties of two molecular-based magnetic helices, composed of 3d metal Co and Mn ions 
bridged by Nitronyl Nitroxide radicals, are investigated by density functional calculations. Their 
peculiar and distinctive magnetic behavior is here elucidated by a thorough description of their 
magnetic, electronic, and anisotropy properties. Metal ions are antiferromagnetically coupled with 
the radicals, leading to a ferrimagnetically ordered ground state. A strong metal-radical exchange 
coupling is found, about 44 meV and 48 meV for Co- and Mn-helices, respectively. The latter have 
also relevant next-nearest-neighbor Mn-Mn antiferromagnetic interactions (of ~ 6 meV). Co-sites 
are characterized by non-collinear uniaxial anisotropies, whereas Mn-sites are rather isotropic. A key 
result pertains to the Co-helix: the microscopic picture resulting from density-functional calculations 
allows us to propose a spin Hamiltonian of increased complexity with respect to the commonly 
employed Ising Hamiltonian, suitable for the study of finite-temperature behavior, and that seems 
to clarify the puzzling scenario of multiple characteristic energy scales observed in experiments. 


I. INTRODUCTION 

Magnetic bistability at the molecular level has at¬ 
tracted broad interest over the last decades in view of 
magnetic-storage applications 1 2. From a fundamental 
perspective, molecules carrying a small number of in- 
teracing paramagnetic centers but behaving like bulk 
magnets, otherwise known as Single-Molecule Magnets 
(SMMs), allowed the direct observation of quantum tun¬ 
neling of the magnetizatiorP, of the Berry phasd- 4 and 
magnetic chiral degrees of freedom 5 (just to cite few 
remarkable phenomena) as well as the manipulation of 
magnetization by light irradiation 6 or via an applied 
voltag^. A large easy-axis magnetic anisotropy is a 
key-requirement to achieve bistability (namely, blocking 
of the magnetization) in molecular magnets and is pro¬ 
vided by metal ions with unpaired electrons carrying a 
finite total angular momentum. These anisotropic build¬ 
ing blocks are coordinated to organic ligands that are 
functionalized in order to favor or shield the propagation 
of exchange interactions. This synthetic strategy nat¬ 
urally leads to the formation of magnetic systems with 
reduced dimensionality^, in which exchange paths may 
have finite connectivity (like in SMMs) or create one- 
(1D) or two-dimensional networks. 

During the 70s and the 80s, molecular ID systems have 
been the object of intense investigation and helped test¬ 
ing the basic princ iples of equilibrium statistical physics 
on realistic systems 10 1=4 The natural evolution of this re¬ 
search line led to study quantum-phase transitions, which 
are still actively explored on ID magnetic systemd^®31. 
Starting from 2001, the observation of slow relaxation 
of the magnetization in molecular spin chains promoted 
them to the role of prototypical systems for the study 
of out-of-equilibrium phenomena as well. By analogy 


with SMMs, slow-relaxing spin chains were called Single- 
Chain Magnets (SCMspP However, the basic properties 
of SCMs are only partially similar to SMMs, since in the 
first ones slow relaxation is not exclusively determined 
by the total magnetic anisotropy. In fact, in SCMs this 
phenomenon may be associated either with the develop¬ 
ment of short-ranged spin-spin correlations upon cooling 
or with the nucleation of a domain wall (DW) at a defect 
sittpHUl 

SCM behavior has been first observed in the 
compound Co(hfac) 2 (NITPhOMe), hereafter called 
CoPhOMe, which is composed of Co(hfac )2 moi¬ 
eties bridged by NITPhOMe radicals (where 
hfac=hexafluoroacetylacetonate and NITPhOMe=4’- 
methoxy-phenyl-4,4,5,5-tetramethylimidazoline-l-oxyl- 
3-oxide) arranged in chiral ID arrays 20 . Then, slow 
relaxation of the magnetization and hysteresis effects 
observed in CoPhOMe were rationalized in terms of 
the kinetic Ising model proposed by Glauber 21 . In the 
comparatively large literature on SCMs that followed, 
slow relaxation has essentially always been interpreted in 
terms of the Glauber model (or variations of it), which 
indeed allows to explain several remarkable features of 
CoPhOMd ^ 20 * 22 *^. Within the Ising model both the 
correlation length and the relaxation time are predicted 
to diverge at low temperature according to an Arrhe¬ 
nius law and with the same energy barrier. However, 
experimentally these two energy scales turn out to be 
different in CoPhOMe 2 -^, thus calling for a deeper micro¬ 
scopic understanding of electronic and thermodynamic 
properties of this compound. In the present study, 
this issue is addressed showing that the mismatching 
between the two energy scales mentioned above can 
be justified relaxing the hypothesis of large (virtually 
infinite) uniaxial magnetic anisotropy, underlying the 
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description in terms of the Ising model. 

According to the present understanding, the emer¬ 
gence of SCM behavior in a given molecular compound 
requires that i) some anisotropy prevents the magnetiza¬ 
tion from reorienting easily, ii) the exchange coupling 
mostly develops along one dimension and Hi) the re¬ 
laxation time of the magnetization becomes macroscopic 
well above the temperature at which residual 3D inter¬ 
actions trigger magnetic ordering. Regarding the first 
requirement, high-spin d 7 cobalt(II) ions in distorted 
octahedral environment are the source of anisotropy in 
CoPhOMe. The ID character (ii) is guaranteed by the 
strong metal-radical antiferromagnetic exchange along 
the chain and by the absence of an exchange path be¬ 
tween magnetic centers belonging to different chain^^S! 
The (residual) dipolar interaction indeed couples also 
spins belonging to different chains of the crystal, but it is 
considered to be too small and partially frustrated to in¬ 
duce 3D magnetic order in CoPhOMe at relatively high 
temperatures. Interestingly, its isostructural Mn-based 
compound, Mn(hfac) 2 (NITPhOMe) (MnPhOMe in the 
following), comprising isotropic high-spin d 5 ions, does 
not show slow relaxation of the magnetization but ex¬ 
hibits a transition to a magnetically ordered phase at 
low temperature (4.8 K) 27 . These isostructural molecu¬ 
lar helices have allowed a comparative investigation of the 
interplay between structural chirality and magnetism, ev¬ 
idencing a giant magneto-chiral dichroism in the hard X- 
ray range for the Co-based noncollinear spin chain, ab¬ 
sent in the isotropic manganese analogue 28 . Moreover 
visible light has been found to promote fast relaxation in 
CoPhOMe through a kick-off mechanism for the nucle- 
ation of domain wallP^. 

A comparative first-principles study of these archetypi¬ 
cal spin chains is therefore of great relevance and may evi¬ 
dence which basic ingredients are needed to produce slow 
dynamics avoiding - at the same time - 3D magnetic or¬ 
dering. As for the last requirement (m), from the micro¬ 
scopic point of view a crucial issue is the relative strength 
of magnetic anisotropy and characteristic exchange en¬ 
ergy, the former being also a key requirement for the ob¬ 
servation of magneto-chiral effects. Here, we address this 
issue by studying the electronic and magnetic properties 
of both CoPhOMe and MnPhOMe via density functional 
calculations. From these, an estimate for the intrachain 
exchange coupling among different magnetic centers is 
provided as well as for the magnetic anisotropy on the 
magnetic ions. The ab—initio analysis, then, allows us to 
propose an effective, albeit realistic, spin model suitable 
to study the thermodynamic properties as well as slow 
relaxation effects of CoPhOMe. Remarkably, our results 
suggest that the exchange interactions in this SCM are 
larger than what previously believed and represent the 
largest energy scale in the model, implying the emergence 
of broad DWs and thus directly affecting the slow relax¬ 
ation properties of the SCMs; moreover, the increased 
complexity of our proposed thermodynamic model can 
qualitatively explain the observation of multiple energy 


scales in CoPhOMe. 

The paper is structured as follows. In the first part 
(Sec. 0, the description of the microscopic properties 
of the systems investigated by density functional calcu¬ 
lations is presented. After a brief methodological sec¬ 
tion repor ting the technical details of the calculations 
(Sec. IIA), the structural properties of the magnetic 
chains are reviewed, with partic ular a ttention given to 
the magnetic building units (Sec. IIB). A detailed anal¬ 
ysis of the electronic properties is given in Sec. |II C| fol¬ 
lowed by a discus sion of the magnetic properties of two 
chains (Sec. IID). The study of the energetics of sev¬ 
eral higher-energy spin configurations allowed to identify 
and estimate the exchange interactions at play, whereas 
the magnetic anisotropy has been quantified in related 
mononucle ar C o- and Mn-complexes. In the second 
part (Sec. Ill), the results obtained by density func¬ 
tional calculations are used as input for an atomistic spin 
Hamiltonian through which the basic thermodynamics of 
CoPhOMe can be described. Finally, the main results are 
summarized in Sec. m 


II. MICROSCOPIC PROPERTIES BY 
FIRST-PRINCIPLES 

A. Methods 

We performed ab-initio calculations within spin- 
density functional theory in the generalized gradient 
approximation 2 ^ as implemented in the VASP materials 
modelling packag^^4 for Co and Mn, the 3d4s states 
were treated as valence states, while for C, N, O, and 
F, the 2s2p were considered as valence states, within the 
projector-augmented waves (PAW) method 32 . Magnetic 
calculations were performed both in the collinear and 
non-collinear formalism!^, including spin-orbit interac¬ 
tion when required. A plane-wave cut-off of 400 eV and 
a 2 x 2 x 1 k-points grid in the primitive cell’s Brillouin 
zone were adopted, and energy convergence was ensured 
by a strict threshold of 10 -6 eV. To model molecular sys¬ 
tems comprising transition metals properly, an accurate 
treatment of electronic correlation is crucial. We have 
addressed this aspect by performing extensive calcula¬ 
tions within the GGA+U approach in the Liechtenstein 
formalism!^, in which a set of two effective parameters 
(U,J) is employed treat electron-correlation, where U de¬ 
scribes the on-site Coulomb repulsion and J is the on-site 
exchange-correction; calculations were also performed in 
the simplified Dudarev’s scheme 3 -^. For each system in¬ 
vestigated, a systematic study of the magnetic properties 
was carried out, spanning U and J parameters over a large 
range of values. In this way, we provided an overview of 
the dependence of the magnetic properties on the corre¬ 
lation parameters, and - most importantly - we ensured 
that our results were robust and conclusions were not 
biased by an inappropriate model description. Unless 
otherwise specified, results for three sets of (U,J) values 
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will be presented in the following, denoted as: set a (Um 
= 5 eV); set b ( U m = 6.5 eV; J m = 1-5 eV) i.e. an ef¬ 
fective U as in set a but with the inclusion of a sizeable J 
value, which turned to be crucial to the aim of modeling 
the magnetic anisotropy; and finally, a set c (U m = 6.5 
eV; J m = 1-5 eV; Uo=Uw = 6 eV) that enhances the 
degree of localization by inclusion of a Hubbard-U also 
on p-states of O and N atoms. 




hfac 


NITPhOMe 



FIG. 1. Color online. Crystal structures of isostructural 
CoPhOMe and MnPhOMe, (a) side and (b) top views. The 
M (Co, Mn) atoms are depicted as large blue spheres, and the 
O, N, C atoms are red, green and brown spheres, respectively. 
For the sake of clarity, F and H atoms are not shown. The 
hfac and NITPhOMe molecules are highlighted in the side 
view; the trigonal symmetry of the chains can be recognized 
from the top view. 


B. Structural properties 

The crystal structures of the two chains in shown 
Fig. 0 composed of alternating M(hfac) 2 moieties and 
NITPhOMe (denoted also with the shortened term NIT, 
in the following) organic radicals arranged in ID arrays 
with helical structures generated by a three-fold screw 
axis. Despite the absence of chiral constituents the com¬ 
pounds form enantiopure crystals, crystallising either in 
the chiral P 3i or P ?>2 space groups. The experimen¬ 
tal structures were used for the calculations, with lattice 
parameters: a = b = 11.294 A, c = 20.570 A for the 
Co-chain; and a = b = 11.281 A, c = 20.846 A for the 
Mn-chain. 

We now discuss in closer detail the structural fea¬ 
tures of the magnetic bricks of the chains, starting with 
CoPhOMe. A view of the molecular cluster is shown in 
Fig- [2] The Co-ion is in a distorted octahedral environ¬ 
ment, coordinated by the O atoms to two hfac ligands 
(03 and 04, and 05 and 06 for the two molecules, see 
Fig. §, while the remaining sites host the terminating 
O atoms of the NITPhOMe radicals in cis position, that 
connect two consecutive Co (hfac) 2 moieties in a chain. 


FIG. 2. Color online. Magnified view of the molecular 
cluster of CoPhOMe, comprising the Co-octahedral complex 
and the NITPhOMe radical bridging two Co-atoms. A similar 
structure characterizes also MnPhOMe. 


The Co-octahedral complex has low symmetry, lacking 
a clear regular pattern. Nevertheless, as an attempt to 
rationalize it, a main tetragonal compression of the oc¬ 
tahedron can be distinguished along the axis connect¬ 
ing the O atoms of the two hfac molecules, the 03 and 
06 sites denoted as apical , whose distance is the short¬ 
est one (4.046 A). This compression indeed corresponds 
to the maximum distortion from the octahedral symme¬ 
try observed in the 03-Co-06 angle, which is 166.62°, 
while the angle formed by the two edging O atoms of 
the cis-coordinated NIT-radicals 01-Co-02 (cfr. Fig. [2| 
is 85.25°. Next, a secondary rhombic-like distortion on 
the equatorial plane of the polyhedron can be revealed 
by the different distances between the vertical O atoms, 
namely 4.147 A for 02-05 and 4.190 A for 01-04. Fo¬ 
cusing on the bridging NITPhOMe radicals, in fact two 
moieties can be distinguished as shown in Fig. [2j de¬ 
noted as a and /3 moieties. These are characterized by 
non-equivalent structural patterns, whereby Co-0 bonds 
are slightly larger in a (Co-Ol=2.108 A) than /3 (Co- 
02=2.097 A) by ~ 0.5%, and viceversa the O-N and 
N-C bonds are shorter in a respect to /? moiety, by 1% 
and 2% respectively. 

Moving to MnPhOMe, the Mn-octahedral complex 
shows a similar distortion pattern to the Co-counterpart, 
characterized by a main tetragonal axial compression, 
with apical O atoms distances of 4.026 A , and a smaller 
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CoPhOMe 



FIG. 3. Color online. Computed density of electronic states 
(DOS) of CoPhOMe (upper panel) and MnPhOMe (lower 
panel). Total DOS are shown (rescaled by a factor 0.1 for 
the sake of clearness) together with the projected DOS rel¬ 
ative to metal (Co, Mn) d-states, and p-states of selected O 
and N atoms of the NITs. 


distortion on the equatorial plane (trans O’s distances 
are 4.326 and 4.252 A ). The NITPhOMe features a sim¬ 
ilar structural pattern to the one discussed above for the 
Co-case, whereby the Mn-0 bonds are larger in the a 
(2.143 A) with respect to the /3 moiety (2.121 A) by 1%, 
whereas the O-N and N-C bonds are shorter in the a 
moiety by 2.5% and 0.1% compared to /3. Globally, Co- 
O bond lengths vary in the range (2.03 - 2.11 A), whereas 
Mn-0 ones are in the range (2.09 - 2.19 A), with Mn-0 
bonds longer than Co-0 ones by about 3 - 4%. 


C. Electronic properties 

Experimentally, CoPhOMe is described by Co(II) ions 
in the high-spin 3/2 state, and NITPhOMe radicals car¬ 
rying a free electron (spin 1/2) that mediates the ex¬ 
change interaction between the Co ions. As already 
stated, Co ions are characterized by a large easy-axis 
magnetic anisotropy 9 20 . Aiming at a systematic analysis 
of the dependence of magnetic properties on the corre¬ 
lation parameters, we initially included the Hubbard-U 
parameter in the d -Co states and varied it in the en¬ 
ergy interval (0 - 5) eV, by searching for the magnetic 
ground state. We found that the inclusion of U Co val¬ 
ues at least of 3 eV is needed to correctly reproduce the 
high spin state of Co(II). As expected, the local magnetic 
moment on Co enhances at higher U values; an optimal 
value of Uco=5 eV was chosen. Let us briefly discuss the 
computed electronic properties of the magnetic chains, 
and their implications for the magnetic properties. In 
Fig. I the total density of electronic states (DOS), to¬ 


gether with the projected DOS of the magnetic elements 
of the chains, the M-ions and the O and N-sites of the rad¬ 
icals, are shown. These calculations refer to the (U,J)-set 
a. For these U values, the chains are small gap insulators, 
with an HOMO-LUMO gap that amounts to 0.60 eV in 
both Mn and Co cases. We however caution the reader 
that the quantitative estimate of the gap value could be 
strongly affected by the exchange-correlation functional, 
as the starting GGA-functional is known to severely un¬ 
derestimate the energy position of the excited states. 
The HOMO and LUMO states are mostly composed of 
the p-states of O and N of the NITs with equal weight, 
while the M-states are pretty far from the gap region. 
In CoPhOMe, the high-spin (S=3/2) electronic config¬ 
uration is here visible in the d-resolved Co-states: the 
spin-up channel is fully occupied, the unoccupied states 
in the spin-down channel are resolved as three peaked 
states, a t^ g state at 3 eV above the Fermi level (E^), 
and two e g states at higher energy (^ 4 eV above Ep). 
The occupied d-states strongly hybridize with the O-p 
states of the NITs in the energy range between 3 and 
1 eV below the Ep. The DOS of MnPhOMe presents 
similar features to the Co-case, with the characteristic 
composition of HOMO and LUMO states by O and N 
atoms of the bridging radicals. The main difference is 
found for the unoccupied down-spin states which appear 
at higher energy, as they start at 4.5 eV above the Ep, 
where two major peaks can be distinguished, composing 
the t 2 g states at 5 eV, followed by the e g states at 6 eV. 
Also in this case the occupied Mn -d states are hybridized 
with the O-p states. 


D. Magnetic properties 

1. Magnetic moments 

The converged magnetic ground state of CoPhOMe has 
a ferrimagnetic (FiM) order, in which the neighboring 
Co-atoms and NITPhOMe molecules are coupled anti- 
ferromagnetically (AF), and well reproduces the picture 
of a Co-ion in high spin state and a free electron delo¬ 
calized on the bridging NITs. A total magnetic moment 
of 6 pb /unit cell is found for this FiM state, as expected 
since the cell can be seen as containing three AF cou¬ 
pled (Co-NIT) magnetic dimers, each one featuring an 
effective moment of 2 pp- In Tab. [I| the values of the 
projected magnetic moments of the Co-ions as well as of 
the NITs are reported. The Co site has a local moment 
of 2.72 pb , and the free electron in the bridging NITs is 
localized on the O-N groups, with a slight unbalance of 
local magnetic moment towards the O-N group in the (3 
configuration (0.04 pp) with respect to the a one; the C 
atom is weakly AF coupled to the two O-N groups. 

Turning to the other chain, in which the Mn(II) ions 
are in 3d 5 electronic configuration, a FiM ground state 
was also found, given by the AF coupling of the mag¬ 
netic bricks Mn and NITs, and a total magnetic moment 
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TABLE I. Projected magnetic moments on M (Co, Mn) and 
on the O-N-C-N-O group of the NITPhOMe radicals corre¬ 
sponding to the FiM magnetic ground state, for CoPhOMe 
and MnPhOMe. The atomic projections are performed within 
spheres of radius as follows: Co (1.30 A), Mn (1.32 A), O (0.82 
A), N (0.74 A), and C (0.86 A). Values are in /is- 


M 

O a 

N a 

C 

N/3 

0/3 

Co (2.72) 

-0.16 

-0.18 

0.05 

-0.20 

-0.17 

Mn (4.59) 

-0.16 

-0.18 

0.05 

-0.18 

-0.16 


amounting to 12 /i^/unit cell. More precisely, the pro¬ 
jected magnetic moments are 4.59 he for the Mn-ions; 
the unbalance between local magnetic moments a ver¬ 
sus (3 on the O-N groups, as discussed above for the Co¬ 
counterpart, here is reduced, with moments of —0.16 and 
—0.18 ii b on the O and N atoms, respectively. 


2. Spin density 

A complementary detailed picture of the magnetic 
structure of the CoPhOMe is provided by the spin den¬ 
sity isosurfaces shown in Fig.[4| with positive spin density 
centered on Co-sites, whereas the negative spin density 
(the radical’s free electron) is equally distributed on the 
p-states of the O and N atoms of the NITs (cfr. Fig.[4|a)). 
The MnPhOMe is characterized by a similar spin density 
distribution (not shown). The only difference that can 
be discerned between the two chains pertains the shape 
of the spin density localized on the M-ions, that mir¬ 
rors their electronic structure. As a matter of fact, the 
Mn-ion are characterized by a perfectly spherical distri¬ 
bution of the spin density, as shown in the 2D sections 
cuts (Fig. Qd)) due to the full (zero) occupancy of the 
spin-up (down) channel. On the other hand, the spin 
density on the Co-sites shows distinctive features of par¬ 
tially occupied e g orbitals, whereby the lobes directed 
towards the O’s of the octahedral complex can be dis¬ 
tinguished in Fig. (4^c). Moreover, the non-collinearity 
of local anisotropies of the Co-octahedral complexes can 
be seen in Fig. J^b), each showing features which are ro¬ 
tated by 120° with respect to the preceding Co-sites by 
trigonal symmetry. 


3. Exchange coupling constants 



FIG. 4. Color online. Isosurface plots of the spin density 
of CoPhOMe: a side view in (a), and top view in (b), the 
latter displaying only the positive spin density for the sake of 
clearness. Positive (negative) values of the spin density are 
depicted as yellow (cyan) lobes. 2D section cuts of the spin 
density taken in a plane parallel to ab and centered on M-ions, 
for CoPhOMe in (c) and MnPhOMe in (d). 
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In order to estimate the magnetic interactions at play 
in the chains, spin-polarized density functional calcula¬ 
tions were performed for an extensive set of collinear 
magnetic configurations. Such configurations, as shown 
schematically in Fig. [5| include higher energy configura¬ 
tions besides the FiM ground state. These configurations 
are obtained by spin flipping of single magnetic units (M- 
and NIT-spin flips), or of dimer M-NIT units in the a and 


FIG. 5. Color online. Schematic view of the set of collinear 
magnetic configurations modeled to calculate the exchange 
interactions in MPhOMe chains, with M = Co, Mn. From left 
to right: the FiM (ground state), the single M- and NIT-spin 
flips, the dimer (M-NIT) spin-flips at the a and f3 bondings, 
and finally the FM magnetic state. The M- and NIT-spins 
are depicted as large blue and small red arrows, respectively. 
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(3 bonding (i.e. frustrating only interactions of the a- or 
/3-type, see Fig. [5|, and include the FM magnetic con¬ 
figuration as well. Note that achieving the convergence 
for these selected magnetic configurations is not a triv¬ 
ial task, as the free-electron on NITs tends to delocalize, 
hence a fine tuning of proper magnetic moments initial¬ 
izations and constraints was needed. The energetics of 
this set of magnetic configurations for CoPhOMe and 
MnPhOMe is reported in Tab. [TI| These results pertain 
to the set a of (U,J) parameters, but consistent values 
have been obtained also for sets b and c (not shown). For 
CoPhOMe, the FM state, with a magnetic moment of 12 
/is/cell, is 572 meV higher in energy respect to the FiM 
magnetic ground state. Consistently, the energy cost of 
single Co- and NIT-spin flips is similar (184 meV). The 
most interesting result concerns the (a or (3) dimer spin 
flips: there is a remarkable energy difference between the 
two configurations, with the a-bond spin flip favored by 
75 meV with respect to the (3 one. Moreover, the (Co- 
NIT) a -spin-flip turns out to be by far the most favorable 
static magnetic excitation. Moving to the Mn-chain, the 
FM state has a magnetic moment of 18 /jlb /cell, and is 
728 meV higher than the FiM ground state. The energy 
cost of a Mn-spin flip is 216 meV, and similarly to the 
Co-case, the (Mn-NIT) a -spin-flip is the favored magnetic 
excitation, while the f3 one is higher in energy by 40 meV. 

An energy-mapping analysis within the broken- 
symmetry approacbP^ was carried out to derive the ex¬ 
change coupling constants, whereby the total energies ob¬ 
tained by first-principles calculations were then mapped 
onto a Heisenberg spin Hamiltonian: H ex = ^m-nit + 
Hnnn comprising nearest-neighbors (M-NIT) and next- 
nearest-neighbors (NNN) contributions: 


^M-NIT 

Hnnn 


$2 r ' S2r+1 + Jp $2r ' S 2r-l] 
r 


[Jm-M $2r • $2r+2 

r 

TJnIT-NIT s 2r-l ' S2r+l] 


(la) 

(lb) 


with S 2 r (even sites) indicating the spins of the metal 
ions and S 2 r +i (odd sites) those of NIT radicals, where 
spins are assumed with unitary modulus in order to di¬ 
rectly compare the strength of the exchange couplings in 
the two chains. The index r labels different metal-NIT 
pairs in each chain. Equation ( [Ta| ) describes the coupling 
between the M-ion and the NIT-radical sublattices, the 
a and /? bondings contributing with a different coupling 
constant, whereas Eq. (lb) takes into account both the 
M-M and NIT-NIT exchange interactions. 

The estimated exchange interactions are shown in 
Tab. Ill Analyzing the CoPhOMe (set a), we first notice 
the large difference in the Co-NIT exchange coupling for 
the two bondings, (3 being stronger by about 19 meV. 
Note also the presence of longer-range NNN interactions 
between Co’s and NIT’s spins of AF-type (of equivalent 
strength, about 1.6 meV). Moreover, by including a finite 
J in the GGA+U Lichtenstein approach on Co-sites (set 


TABLE II. Computed relative energies of the collinear mag¬ 
netic configurations for CoPhOMe and MnPhOMe, compris¬ 
ing the FiM, the M- and NIT-spin flips, the dimer (M-NIT)- 
spin flips, and the FM. Calculations performed with Um=5 
eV (set a), energy values are in meV. 


Magnetic config. 

CoPhOMe 

MnPhOMe 

FiM 

0.0 

0.0 

M-spin-flip 

184.1 

216.0 

NIT-spin-flip 

184.5 

216.6 

(M-NIT) a -spin-flip 

139.3 

173.4 

(M-NIT) ^-spin-flip 

214.4 

213.7 

FM 

572.2 

727.7 


TABLE III. Exchange-coupling constants in CoPhOMe and 
MnPhOMe, as derived by Eq. 0 and with the convention of 
normalized spins in magnitude, computed for the three sets 
of (U,J) correlation parameters. Values are in meV. 


CoPhOMe 

J OL 


J Co-Co 

Jnit-nit 

a 

-34.8 

-53.6 

-1.7 

-1.6 

b 

-31.3 

-45.9 

-1.3 

-1.4 

c 

-23.8 

-34.7 

-0.8 

-0.4 

MnPhOMe 

J OL 

J0 

^Mil-Mil 

Jnit-nit 

a 

-43.3 

-53.4 

-6.6 

-6.5 

b 

-43.7 

-53.9 

-6.7 

-6.4 

c 

-36.5 

-44.3 

-5.0 

-4.3 


5), and by increasing further the correlation on O and N 
atoms (set c), the exchange-coupling constants are sub¬ 
ject to a reduction, as expected. We note that a uniform 
reduction for the NN-interactions is found, expressed by 
a constant ratio J a /Jp ^0.7 for the three (U,J)-sets here 
considered. However, this is not the case for the NNN- 
interactions, whereby the Jnit-nit is more significantly 
reduced by the inclusion of U on the O- and N- atoms 
than the Jco-Co- This can be understood considering 
that the exchange coupling Jnit-nit involves two NITs 
and an intermediate Co-atom, and it is therefore more af¬ 
fected by an increased localization onto the NITs atoms. 
The value of the Co-NIT exchange coupling constant here 
calculated is larger than the experimental estimation ob¬ 
tained by magnetic susceptibility measurements 20 . How¬ 
ever, it is important to emphasize that the experimen¬ 
tal estimates of the exchange couplings strongly depend 
on the spin Hamiltonian that is assumed to compute the 
magnetic susceptibility. This is crucial for CoPhOMe, for 
which the assumptions of a pure Ising model turns out 
to be inadequate, as it will be shown in Sec. IH where 
the thermodynamic properties of the system are derived 
by using a more sophisticated effective spin-Hamiltonian 
based on density functional results. Note moreover that 
these exchange energy values are consistent with ab initio 
calculations within Quantum Chemistry approach car- 
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ried out on a simplified cluster model composed of a rad¬ 
ical Co-NIT pair (in the a-bonding) 26 . The authors ob¬ 
tained a value J a -163.5 cm 1 (-20.3 meV), evaluated 
from the energy gap between the two lower spin-states, 
i.e. the triplet (S=l) and the quintet (S=2) spin-states 
corresponding to the AF and FM coupling of the Co-NIT 
spin-pairs, respectively. 

Moving to MnPhOMe, the computed exchange inter¬ 
actions are AF-type, both Mn-NIT and NNN ones, and 
are considerably stronger than those found in the Co¬ 
chain. Particularly, the a bondings are stronger than 
the corresponding Co ones by about 10 meV. Impor¬ 
tantly, the NNN interactions are here sizable (about 6.5 
meV). Similarly to CoPhOMe, the magnetic exchange at 
the a bond is weaker than at the /3 one, with a consis¬ 
tent ratio Ja/J/3 ~ 0.8 for all the considered sets. Ex¬ 
perimentally, among the Mn-NIT chain complexes, Mn¬ 
PhOMe has been characterized as the compound with the 
largest value of antiferromagnetic coupling constant 27 . 
By using a NN Heisenberg model to reproduce the ex¬ 
perimental susceptibility, an exchange coupling constant 
^Mn-NIT = —344 cm -1 (—42.6 meV) was estimated. Due 
to the negligible anisotropy of Mn(II)-ions, the g factor of 
both Mn and NIT spins was fixed to 2.0 value for the fit¬ 
ting. Our results are in good agreement with this exper¬ 
imental value, and fully justify the choice of an isotropic 
Heisenberg model (see below). 

Some considerations about the calculated exchange 
coupling constants are here given, (i) The magnetic 
building units in both compounds are chemically rather 
similar, i.e. Co or Mn ions coupled via direct type mag¬ 
netic exchange to nitronyl-nitroxide radicals; the main 
difference consisting in an enhanced exchange pathway 
for the Mn-compound, due to the half-filled d-shell elec¬ 
tronic configuration. Therefore, our results of exchange 
coupling constants of strength of the same order of mag¬ 
nitude for both complexes, yet with a stronger coupling 
for the MnPhOMe, are fully consistent based on these 
qualitative arguments, (ii) The strength of the exchange 
coupling strongly depends on the structural parameters 
of the magnetic constituents. As a matter of fact, a 
change in structural parameters such as bond lengths 
and angles will modify the degree of overlap of the lo¬ 
cal magnetic orbitals contributing to the exchange path¬ 
ways, and this will impact on the resulting strength of 
the exchange coupling. Recently, the magneto-structural 
correlations between Mn(II) ion and NITs in trans- and 
czs-coordinated model complexes were investigated by 
broken-symmetry DFT approactl 37j . The authors found 
that the exchange coupling constant Jmii-nit increases 
linearly by decreasing the bond distance Mn-0 (O of 
the NIT) in these complexes. If we consider solely the 
effect of the Mn-0 bond length variation going from 
<a-bond (2.14 A) to the (3 -one (2.12 A) in MnPhOMe, 
from their results we can roughly estimate an increase 
A JMn-NiT of ~ 6%. Other structural parameters are in¬ 
volved as well, and that may account for the variation 
of ~ 20% for the two bondings in MnPhOMe, as indi¬ 


cated by our calculations. In this regard, bond angles 
variations are even more important, as they can induce 
variations in the exchange coupling constants up to 40 
- 509P 1 . (in) Finally, it is worth emphasizing that the 
method adopted - by considering the whole experimen¬ 
tal structure within a periodic boundary approach and 
examining various magnetic configurations - allowed us 
to enrich the picture of the magnetic interactions at play: 
quantities which are difficult to access experimentally 
and neglected within a cluster-model approach. Notably, 
although metal-ions are far apart (the intra-chain Co-Co 
and Mn-Mn distances are 7.815 A and 7.878 A, respec¬ 
tively), we do appreciate non-negligible NNN-exchange 
interactions, which overall may be particularly relevant 
for the specific case of MnPhOMe. Interestingly, in anal¬ 
ogous chains comprising lanthanide ions and NIT radicals 
NNN interactions are dominating^. 


4 • Magnetic anisotropy 

The anisotropy of the Co-compound is ultimately of 
magneto-crystalline type, namely single-ion anisotropy 
deriving from the dJ electronic configuration in the oc¬ 
tahedral field, with a hole in the t^g orbitals. This elec¬ 
tronic structure yields a sizable orbital magnetic moment 
on the Co-sites, and this leads to strong anisotropy due 
to the large first-order spin-orbit coupling. The ion Co 2+ 
is therefore susceptible to the anisotropy of the chem¬ 
ical environment, whereas this is not the case for the 
Mn 2+ , which is in d 5 configuration and with negligible 
orbital angular momentum. This leads to the different 
magnetic behavior of the two systems. Due to the low 
symmetry of the metal-sites local chemical environments 
and their distortions from octahedral symmetry, a direct 
calculation of the anisotropy for the real chemical sys¬ 
tems is needed. The characterization of the local mag¬ 
netic anisotropy of the M-sites in the chains turns out 
to be unfeasible due to the impossibility to single out 
the single-ion anisotropy from the exchange interactions 
with the confining radicals. To this purpose, we have 
hence resorted to the molecular compound Co(hfac) 2 - 
(NITPhOMe) 2 , hereafter called C 0 NIT 2 , which is essen¬ 
tially the mononuclear variant of the chain, containing 
the octahedrally coordinated Co(II) ions with two hfac 
molecules and two NITPhOMe ligands. The latter are 
coordinated in cis configuration to the central Co-ion 
through a single ON groupP^. 

Since the metal ion has essentially the same local co¬ 
ordination environment in the mononuclear compound 
and in the chain, the two complexes are expected to be 
characterized by the same local anisotropy properties. 
The Co single-ion anisotropy in the compound C 0 NIT 2 
was evaluated by performing magnetic calculations (in 
the non-collinear formalism and inclusive of SOC) for 
the spin-trimer magnetic structure, i.e. with the spins of 
the NITs collinear among them and AF-coupled to the 
spin of the intermediate Co-atom. In this way, the two 


Co-NIT exchange interactions are fully satisfied and the 
energy variations obtained by rotating the orientation 
of the (total) spin-trimer axis can be safely attributed 
entirely to the Co single-ion anisotropy, since the spins 
of the attached radicals are isotropic. The experimen¬ 
tal structure of the C 0 NIT 2 , which crystallizes in the 
triclinic PI space group^, was considered for the calcu¬ 
lations. The Magnetic Anisotropy Energy (MAE) was 
characterized by total energy calculations with the axis 
of the total spin-trimer spanning the three planes of a 
Cartesian system XYZ (in steps of 15°). 

Regarding the Mn-compound, the experimental struc¬ 
ture is not available in literature. Therefore, we used 
the structure of the C 0 NIT 2 and replaced the central Co 
atom by Mn. The Mn-0 bond lengths in the octahe¬ 
dral complex are shorter in this model system compared 
to the experimental structure of the MnNIT 2 compound, 
and this is expected to slightly overestimate the magnetic 
exchange between the Mn-ion and radicals. However, be¬ 
ing aware of this limitation of the model, we resort here 
to the molecular cluster solely to inspect the single-ion 
anisotropy, and we believe this should not affect signifi¬ 
cantly the results for this specific purpose. 

As an additional methodological remark, we found that 
the magnetic anisotropy critically depends on the pa¬ 
rameter J in the GGA+U formalism, while the role of 
the parameter U is negligible. Specifically, whereas the 
anisotropy profile is overall consistent, with location of 
minima and maxima in the energy landscape being un¬ 
altered, the strength of the anisotropy enhances with in¬ 
creasing J values. Based on that, the results here pre¬ 
sented are relative to set c, with M-ions subject to a 
sizable J value of 1.5 eV. As a matter of fact, calcula¬ 
tions with J=0, i.e. set a, show a reduced anisotropy 
profile. It is worth noting that non-collinear calculations 
were performed by achieving full self-consistency after 
the inclusion of SOC-term with a more accurate energy 
threshold, as we realized that standard magnetic calcula¬ 
tions with non-selfconsistent SOC cycle yield anisotropy 
properties largely quenched. The non-collinear magnetic 
configurations were generated by adopting the constrain¬ 
ing technique of the local magnetic moments orientation, 
as implemented in the VASP code. This approach is es¬ 
sentially based on the inclusion of an energy penalty term 
to the total energy that drives the local moments to the 
desired directions. We set the energy threshold to 10 -7 
eV and, after proper initialization of the magnetic mo¬ 
ments along defined spatial directions of the total spin, 
increased stepwise the weight of the penalty term, until 
convergence of the non-collinear constrained SOC calcu¬ 
lations was achieved. 

MAE profiles for both compounds are shown in Fig. [6] 
It is evident that, while the Mn-compound is nearly 
isotropic, with the largest energy variation of about 0.15 
meV, the CoNIT 2 shows a sizable magnetic anisotropy 
with the largest energy barrier ^ 3.85 meV (46K). Addi¬ 
tionally, upon closer inspection, the anisotropy on each 
Co-ion turns out to be uniaxial. This is evidenced 


by the markedly deeper energy barrier for magnetiza¬ 
tion reversal experienced by the spin-trimer in the XY- 
plane compared to the shallower barriers (saddle points) 
along the other two sampled planes. These results are 
in good agreement with the experimental characteriza¬ 
tion of the MnPhOMe as an isotropic Heisenberg-like 
FiM spin chain, as well as with the description of the 
CoPhOMe as a non-collinear FiM spin helix, deriving 
from uniaxial character of Co-octahedral complexes. In 
this second spin chain, the easy axis of each Co is not 
collinear with respect to the ones of neighboring Co¬ 
atoms but follows the three-fold periodicity of the crystal. 
By mapping the estimated local easy-axis of the Co ion 
in the molecular compound onto the crystal frame of the 
chain, the tilt angle formed by the Co spins relative to 
the chain axis c (polar angle 0) can be estimated, yield¬ 
ing 0 = 50 - 55°. This value is in reasonable agreement 
with the experimental estimate given in Ref. [40] (see be¬ 
low). However, a direct comparison of the MAE obtained 
by DFT-calculations with the EPR-characterization of 
the anisotropy 39 in the compound C 0 NIT 2 shows some 
inconsistencies regarding the absolute orientation of the 
anisotropy easy-axis. This aspect requires further inves¬ 
tigation and will be the subject of future work. Notwith¬ 
standing, it is worth emphasizing that this discrepancy 
between theory and experiments about the absolute easy- 
axis direction in the mononuclear compound will impact 
- once mapped onto the chain crystal frame - only on the 
azimuthal angle (i.e. the orientation on the ab plane), 
which is overall not relevant due to the chain trigonal 
symmetry, while the polar angle 0 is well reproduced. 

Additionally, the anisotropy was assessed directly in 
the crystal frame of the chain compound. A model sys¬ 
tem was used for this purpose in which two Co atoms 
were replaced by Zn’s in the unit cell of CoPhOMe. In 
this way, the diamagnetic Zn(II) ions turn off the mag¬ 
netic interactions with the neighboring radicals, and the 
magnetic complex is hence simplified to a spin-trimer 
comprising the Co-site AF-coupled to two lateral NITs, 
though in the helical structure. By rotating the axis of 
the total spin of the complex along defined directions, the 
anisotropy of the local Co-octahedral environment was 
probed. The set of directions included the Co-0 axes 
and bisectors of the octahedral local frame. Due to the 
demanding computational task required, the SOC term 
was included in a non-self-consistent way. Although the 
comparison remains qualitative since non-self-consistent 
calculations do not account for the whole anisotropy 
strength, the results are in accord with those performed 
on the mononuclear-complex, namely easy-axis as well as 
energy valleys/rises of the MAE profile in the two local 
frames coincide. 

Static measurements of the magnetization in 
CoPhOMe show that the magnetization relaxes slowly 
below a temperature of 6 K when an external magnetic 
field is applied parallel to the trigonal axis, while no 
(dynamic) hysteresis is observed for magnetic fields 
applied on the trigonal plane 20 . This observation can be 
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Plane YZ 



FIG. 6. Magnetic Anisotropy Energy profiles for Co- and 
Mn-based molecular compounds C 0 NIT 2 and MnNIT 2 , as the 
orientation of the total spin is varied along three planes of 
a Cartesian system XYZ. The axes (X,Y,Z) of the reference 
Cartesian system are related to the crystallographic directions 
(a, b , c) of the triclinic lattice as follows: a = X; b lies on the 
XY-plane at 64.51° with respect to X-axis, the c-axis is tilted 
of 11.97 0 with respect to the Z-axis, and with projection 
on the XY-plane lying on the first quadrant. Energies are 
relative to the value on the X-axis direction. 


ascribed to non-collinearity among anisotropy axes, as it 
will be discussed in Sec. mm 


III. FINITE TEMPERATURE PROPERTIES 

The thermodynamics of spin chains is dictated by the 
temperature dependence of the correlation length £, i.e., 
the characteristic scale of decay of pair-spin correlations. 
For infinite chains, this intrinsic length scale is propor¬ 
tional to the product of the static susceptibility (x) mea¬ 
sured in zero magnetic field and temperature: ^ £• 

The last relation is usually employed to extract infor¬ 
mation about the spin-Hamiltonian parameters directly 
from experimental susceptibility data. In the absence 
of anisotropy (as in the case of MnPhOMe) the corre¬ 
lation length is expected to diverge like J ex /T with de¬ 
creasing temperature (J ex being some effective exchange 
constant describing the coupling between nearest neigh¬ 
bors). In spin chains with uniaxial anisotropy, one ex¬ 
pects £ ^ e A ^/ T , where the characteristic energy scale 
may have a residual dependence on T arising from 
renormalization due to spin waved^ (&b = 1 will be as¬ 
sumed henceforth). Such a mechanism is not effective 
when elementary excitations are sharp (Ising) DW& 4 ^, 
which happens when the uniaxial anisotropy is compa¬ 


rable to the exchange interaction or larger. In this case, 
= 2 J ex throughout the whole range of temperatures 
where £ exceeds some lattice units (assuming a unitary 
spin modulus). When the exchange energy is larger than 
the anisotropy energy, elementary excitations are broad 
DWs and the relationship between A^ - accessible in 
experiments - and spin-Hamiltonian parameters is not 
straightforward anymore. The anisotropy energy com¬ 
puted for the CoNIT 2 and the estimates of J a and Jp 
obtained for CoPhOMe suggest that in this spin chain 
elementary excitations should actually consist of broad 
DWs. The scenario is further complicated in CoPhOMe 
by the non-collinearity among local anisotropy axes. In 
the following, the question of how the energy barrier A^ 
is affected by the degree of non-collinearity will be ad¬ 
dressed. 


A. Model and results 

To model the thermodynamics of CoPhOMe we pro¬ 
pose the Hamiltonian 

7~L = ^ ^ [^ex ' (&p,2r-\-l A &p,2r— l) 

p,r 

+ Sp,2rD2 r Sp, 2r + I^B (#B • <T p ,2r+l + G B Sp ? 2r) 

( 2 ) 

where Co spin operators are replaced by classical vec¬ 
tors | Sp ? 2 r | = 1 and <7 p ,2r+i are Pauli quantum oper¬ 
ators, representing NIT spins. A uniaxial anisotropy 
was assumed for each Co in the corresponding local 
frame, namely D x = D y = 0 and D z = 46 K (ob¬ 
tained from DFT calculations on C 0 NIT 2 ). Each of these 
anisotropy tensors is reported to the crystal frame by 
a different rotation parameterized by a standard ma¬ 
trix 7£(<£, 0,-0), where cp,6 ,-0 are Euler angles. More 
explicitly D 2r = 0, ip) dia,g(D x , D y , D z )H T (<p, 9, ip), 

where diag (D x , D y , D z ) indicates the 3x3 diagonal ma¬ 
trix with eigenvalues D x , D y and D z . The crystal sym¬ 
metry imposes ip = 2ttv/3 with r = 1,2,3 so that the 
index r labels the atoms inside the cell, p is the cell in¬ 
dex and ideally extends to infinity. Due to our choice 
of vanishing D x and D y , the tensor D^ r is independent 
of the 'ip angle. Therefore, the degree of non-collinearity 
is uniquely parameterized by 6: For 0 = 0° anisotropies 
are collinear and point along the chain axis c; while for 
0 = 90° local anisotropy axes lie on the ab plane, with 
neighboring axes forming an angle of 120°. The coupling 
between Co ions and the magnetic field should properly 
be described by a Lande tensor. However, since the Zee- 
man term does not affect the energy barrier A^, a scalar 
Lande factor was assumed for both Co atoms and NITs, 
with G = 2 x 3/2 = 3 and g = 2 x 1/2 = 1 (including 
the spin modulus). Moreover, only the exchange cou¬ 
pling between each Co and the neighboring radicals was 
considered and set to J ex = —500 K. This value is in¬ 
termediate between J a = —404 K and Jp = —622 K, 
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obtained from the set a of DFT calculations. In Ap¬ 
pendix [A] it is shown that distinguishing between J a and 
J$ does not alter the thermodynamic properties of this 
model in the temperature range of interest and for this 
choice of parameters. The mapping of the model with 
J a ~/~ Jp into a model with a unique exchange coupling 
J ex relies on assuming broad DWs. Therefore, for values 
of D z consistent with the formation of sharp DWs this 
equivalence is not guaranteed. 

Equilibrium properties at finite temperature can be de¬ 
duced from the knowledge of the partition function Z, 
which is related to the free energy F by the fundamental 
equation F = — T In (Z). The magnetization and the sus¬ 
ceptibility are given by the first and the second derivative 
of F with respect to the applied field. These quantities 
have been computed for the infinite chain, i.e. N p oc, 
with the transfer-matrix technique. Details about this 
method, including the treatment of quantum operators 
0>,2r+i 5 are given in Appendix |B| 



FIG. 7. Color online. Susceptibility along the chain axis 
X|| = dM\\/dB\\ computed in B = 0 for 6 = 0° (dots), 40° 
(crosses), 50° (stars), 52° (empty squares), 54° (full squares). 
The arrow indicates the direction of increasing 6. Solid lines 
represent the linear fittings of low- and high-temperature bar¬ 
riers: A^ ow (black line) and A^ lgh (red line). 

We focused on the susceptibility X|| = dM\\/dB\\ along 
the chain axis c computed in zero field for different val¬ 
ues of the angle 0. In Fig. [7] these results are plotted 
as X|| T versus 1/T in a log-linear scale in order to high¬ 
light the exponential divergence of the correlation length 
~ e A ^/ T and possible deviations from it. For each value 
of 0, two “linear regimes” - corresponding to two val¬ 
ues of A^ - can be identified at low and intermediate T. 
For a given 0, the value of A^ is systematically larger at 
low T than at hight T. Black and red lines are drawn 
on top of the X|| T curve computed for 0 = 54° to ex¬ 
emplify the two linear fittings. Enlarging the horizontal 
scale two distinct slopes can be highlighted for smaller 


6 angles as well. The two energy barriers correspond¬ 
ing to each 0 are plotted with symbols in Fig. [8j blue 
stars for the low-temperature Aj? w and red crosses for 

the high-temperature A^ lgh . As non-collinearity among 
anisotropy axes is increased, by increasing 0, both high- 
and low-temperature values of A^ decrease: the former 
passes from 210.4 K in the collinear case (0 = 0°) to 30.9 
K for 0 = 54°; the latter ranges from 280.2 K to 50.9 K 
when 0 is varied in the same interval. The existence of 
two energy barriers Aj? w and A^ lgh is not due to non- 
collinearity. In fact, this phenomenon occurs also in fer¬ 
romagnetic Heisenberg chains with collinear anisotropy 
axes for ratios of D z to the exch ange en ergy compatible 
with the formation of broad DWP 41 42 . As anticipated, 
this is in contrast to what happens in the limit of sharp 
DWs, where A^ takes a single value, right equal to the 
DW energy. Moreover, in Heisenberg chains with broad 
DWs and collinear anisotropy A^ is always smaller than 
the DW energy. At low temperature this is due to the 
interplay between DW excitations and spin waves. More 
precisely, for each DW that is added to the system, two 
Goldstone modes appear in the spin-wave spectrum— 
that are associated with translational invariance of the 
DW center and with the degeneracy with respect to the 
azimuthal angle (e.em Bloch and Neel DW have the same 
energy if D x = D y Besides this, the presence of DWs 
also modifies the density of states of spin waves with 
respect to the case in which the last ones are superim¬ 
posed to a uniform spin profile. All this results in an en- 
tropic contribution that affects the correlation length and 
makes the low-temperature A^ ow smaller than the DW 
energjPHS Still in collinear chains, the further suppres¬ 
sion of A^ occurring at intermediate temperatures has 
been justified in the framework of P olyako v renormaliza¬ 
tion of spin Hamiltonian parameter d 41 * 42 *. 

By analogy with the collinear case, we may infer that 
the suppression of A^ with increasing temperature ob¬ 
served in our non-collinear model - at fixed 0 - also 
emerges from the interplay between spin-wave and DW 
excitations. In fact, with some additional hypotheses 
based on J ex D z (see Appendix [A]) and setting B = 0, 
Hamiltonian ([ 5 ]) can be mapped into the following func¬ 
tional 



-D z (3cos 2 9 — l) ( S Z Y 


dz . 


(3) 

This functional is formally the same as the continuum 
version of the Heisenberg Hamiltonian with collinear 
anisotropy axes discussed before. Note that the penal¬ 
ization due to the misalignment of neighboring spins 
is halved with respect to the ferromagnetic Heisenberg 
chain (|J ex |/4 instead of J ex /2), because the coupling 
between two consecutive classical (Co) spins is medi¬ 
ated by the quantum spin of the interposed NIT, and 
the uniaxial anisotropy is replaced by the 0-dependent 
term D z (3 cos 2 0 — l) /2. Thanks to this mapping, the 
analytic expression giving the cost to create a DW in 
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the Heisenberg chain with collinear anisotropies 46 can be 
used to get £dw(9) = \/2\Jex\D z (3 cos 2 6 — 1). Remark¬ 
ably, the effective anisotropy (and £dw{0) with that) van¬ 
ishes at the magic angle 0 = 54.7°, above which both the 
magnetization and the susceptibility become larger on 
the ab plane than along the c axis (parallel to z). £dw($) 
is plotted as a solid black line in Fig. [8] In collinear 
chains both high- and low-temperature are smaller 
than £dw but proportional to ifM According to Fig. [8j 
this trend seems to be maintained in CoPhOMe, in spite 
of its non-collinear structure. To better verify this trend 
we rescaled £dw(9) to the high- and low-temperature val¬ 
ues taken by A^ in the collinear case (0 = 0°). The 
corresponding curves are plotted with lines of the same 
colors as the fitted barriers A^ (symbols). The agreement 
is quantitative within 5% up to 0 = 45°, while larger de¬ 
viations appear when 0 approaches the magic angle. For 
6 < 45° we can, therefore, express A^ in terms of the 
spin-Hamiltonian parameters and the 9 angle as follows: 

A“ gh ~ 0.69£ dw (6>) = 0.69v/2|J ex |D 2 (3cos 2 6»-l) 

? ;_ (4) 

A\ ow ~ 0.92£ dw (6>) = 0.92 v / 2|J ex |L> z (3cos 2 6»-l) 

where the superscripts stand for high- and low-T, respec¬ 
tively. 



FIG. 8. Color online. Energy barriers A^ ow (blue stars) and 
A^igh ^ re( j crosses) extracted from the susceptibility curves at 
low and intermediate temperatures, as described in the text, 
are plotted with symbols as a function of 0. Error bars are de¬ 
termined by the uncertainty on the fitted slopes. The upper, 
black line represents the DW energy £dw($); the dashed blue 
and the red lines represent the same quantity rescaled to the 
value of A^ obtained for 0 = 0 at low and high temperature, 
respectively (see Eq. @). 

For the sake of completeness, we mention that Hamil¬ 
tonian (2| is suitable to model MnPhOMe as well, pro¬ 
vided that the anisotropy term is removed. In zero 
applied field the relevant observables can be computed 
analytically 47 . 


B. Comparison with experiments 

The angle that anisotropy axes form with the c axis 
in CoPhOMe is estimated to be 0 ~ 50° or larger by 
DFT calculations. A similar estimate was already re¬ 
ported in the early works on this compound based on the 
(small) difference between the saturation magneti zatio n 
measured along the chain axis and on the ab plan d 4Q | 4 ^ . 
For this value of 6 and the computed parameters (set a) 
the DW energy is £dw(50°) = 105 K. In the following we 
will focus on the high-temperature A^ lgh , which is the 
only barrier accessible experimentally. In fact, finite-size 
effects - jnduced by the presence of naturally occurring 
defectd^ 7 ^ - and slow dynamics prevent the correlation 
length, and \ consequently, from diverging indefinitely as 
temperature is lowered. The numerical transfer-matrix 
calculation yields A^ lgh = 76 K, which falls in the cor¬ 
rect range, even if a more quantitative comparison of this 
thermodynamic model with experiments would require a 
direct fitting of the susceptibility data. At present, it is 
important to stress that the observation of A^ of the or¬ 
der of hundred K does not necessarily imply that | J ex | be 
of the order of 50 K, as the Ising model would prescribe: 
With the more realistic model Hamiltonian in Eq. © an 
estimate of J ex = —500 K obtained from ab-initio calcu¬ 
lations can be absolutely consistent with a much smaller 

When it was firstly synthesized, the novelty of 
CoPhOMe consisted in displaying slow relaxation of the 
magnetization. The relaxation time was observed to 
obey an Arrhenius law r = roe Ar / T over ten decades, 
which was expl ained in the framework of the Glauber 
modeP^ 19 1 23 E 25 ! Strictly speaking, the Glauber model 
is a kinetic version of the Ising model, which assumes a 
very large (ideally infinite) anisotropy. Even if in realistic 
SCMs the anisotropy is finite, this description has proven 
to be appropriate at the condition that elementary exci¬ 
tations are sharp DWs. In particular, at low temperature 
the relaxation process initiates by reversing a spin lying 
at one of the two edges of each open chain (see below), 
with an activation energy equal to the DW energy and, 
therefore, to A^. The diffusion of a sharp DW requires 
an addit ional activation energy (usually named A^) at 
each stepP- 6 ^ . Finally, in the large-anisotropy limit, one 
has A r = A^ + A^ 4 . 

The fact that DFT calculations yield |J ex | D z calls 
for a revision of the mechanism behind slow relaxation 
and the relative energy scales, because if elementary exci¬ 
tations are broad DWs the Ising and the Glauber models 
are not applicable. So far there is less quantitative agree¬ 
ment between theory and experiments on SCMs with 
broad DWs and the interpretation of available data is 
still controversial. However, theoretical arguments sug¬ 
gest that the diffusion of a broad DW does not proceed 
by thermal activatiorP^ (i.e. A^ = 0). The energy 
cost to nucleate a DW from one edge of an open spin 
chain is expected to be of the order of the DW energy at 
zero temperature £dw(9). A later experimental study on 
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CoPhOMe - in which the amount of non-magnetic impu¬ 
rities was varied in a controlled way - confirmed that the 
relaxation precess indeed initiates at one edge of a finite 
segment of coupled spins 1719 . Since in CoPhOMe 0 is 
certainly larger than 45°, Eq.Q does not hold exactly. 
However, this equation can still help estimate the rela¬ 
tion between A^ lgh and A r , the latter beeing identified 
here with the DW energy: one expects the measured A^ 
to be roughly 70% of A r . From the value of J ex used 
in Ref. 20] to fit the susceptibility data, we deduce that 
A^ 110 K for that specific sample; in the same refer¬ 
ence it is reported that A r = 154 K, which yields a ratio 
between these two energy barriers in fair agreement with 
our expectation. In conclusion, this alternative paradigm 
based on broad DWs seems plausible. Further details 
about how the diffusion of a broad DW can be mapped 
into the Glauber model are given in Appendix [Cj 

The Glauber model successfully justifies the obser¬ 
vation of slow relaxation in CoPhOMe only for mag¬ 
netic fields applied along the chain axes. As already 
mentioned, this originates from non-collinearity among 
anisotropy axes 24 40 . By virtue of Hamiltonian ©, 
the same phenomenon can be justified in terms of the 
behavior expected for Heisenberg chains with collinear 
anisotropy. As far as 0 is smaller than the magic an¬ 
gle, the effective anisotropy D z (3 cos 2 6 — l) /2 is posi¬ 
tive and the global easy axis points along c. The creation 
of a DW by thermal activation is thus required to reverse 
the magnetization along the c axis (|| z). Conversely, ro¬ 
tating the magnetization on the ab plane should not re¬ 
quire any activation energy. For larger 9 , the spin chain 
becomes easy-plane, but still the magnetization rever¬ 
sal on the ab plane should not be a thermally activated 
process because in this case the system should behave 
similarly to the XY model. 


IV. CONCLUSIONS AND PERSPECTIVES 

In this work, we have addressed - by means of density 
functional calculations - the study of two archetypical 
organic magnetic chains, composed of M(hfac )2 moieties 
bridged by NITPhOMe radicals, with M=Co, Mn. They 
are structurally very similar, but their magnetic behavior 
is remarkably different. CoPhOMe chain is characterized 
by SCM behavior, whereas MnPhOMe is described as 
an isotropic spin chain. A detailed investigation of the 
magnetic, electronic, and anisotropy properties has been 
performed, enriching the understanding of these systems, 
with particular emphasis on the description at the atomic 
scale. 

Consistently with experiments, we find that both 
chains have a ferrimagnetic ground state, as a result of 
the local AF-coupling between the metal-ion and the rad¬ 
ical. The exchange interaction M-NIT is rather strong, 
of the order of 30 - 50 meV for CoPhOMe, and even 
of greater strength for MnPhOMe. The latter is char¬ 
acterized also by sizable next-nearest neighbor exchange 


interactions (between Mn-Mn and NIT-NIT couples) of 
AF-type. Importantly, we find that there are two types 
of M-NIT exchange couplings, of considerably different 
strength. This finds a correspondence in their structural 
features, whereby the bridging NITs are characterized 
by two non-equivalent bondings with the confining metal 
ions. The difference in CoPhOMe is about 19 meV, while 
in MnPhOMe it is ~ 10 meV . 

The electronic properties are pretty similar in the two 
chains, both systems behaving as small gap semicon¬ 
ductors. The HOMO and LUMO states are essentially 
composed of the p-states of the spin-polarized O and N 
atoms of the NITs, while the occupied d-states of the 
transition metals are quite far from the gap, strongly hy¬ 
bridizing with the O states. The magnetic anisotropy 
has been characterized in the mononuclear variant of 
the chains, the M-NIT 2 compound, which is a magnetic 
spin-trimer. We find that the Mn-compound is almost 
isotropic, while the Co-compound is characterized by a 
local uniaxial anisotropy with an estimated energy bar¬ 
rier of the order of 4 meV. This scenario is consistent 
with the phenomenological description of CoPhOMe as a 
non-collinear spin-spiral, with local anisotropy easy axes 
tilted with respect to the helix axis. 

We have proposed a spin Hamiltonian, Eq © , that en¬ 
ables the comparison of magnetic properties determined 
by ab-initio calculations with experiments in a very effi¬ 
cient way, still keeping the rich complexity of CoPhOMe 
(the same Hamiltonian without the anisotropy term de¬ 
scribes also the thermodynamics of MnPhOMe). We 
have focused on the dependence of a characteristic energy 
scale A^ on the structural properties of CoPhOMe, pa¬ 
rameterized by the angle 0 formed by the local anisotropy 
axes of Co ions with the c axis. We derived an analytic 
formula that describes the non-trivial dependence of A^ 
on the MAE, on the Co-NIT exchange coupling and on 
0. This proves the increased complexity of our thermody¬ 
namic model with respect to the Ising model - commonly 
employed to rationalize SCMs - in which it is trivially 
2Jex- 

The relative strength of the Co-NIT exchange coupling 
with respect to the MAE is consistent with the formation 
of broad DWs as elementary excitations. In this sense, 
DFT results provide a new framework to describe the 
static and dynamic properties of CoPhOMe and solve the 
apparent contradiction of A^ being significantly different 
from A r (barrier of the relaxation time): the latter is 
expected to be of the order of the DW energy £dw(#) and 
the former roughly 70% of it. 

From a methodological perspective, the present study 
shows how the molecular approach spontaneously leads 
to isolate, from an extended system, individual build¬ 
ing blocks that can be studied with deeper accuracy, 
both theoretically and experimentally. The crucial in¬ 
formation extracted from those building blocks and from 
their arrangement in the extended system is retained at a 
coarse-grained level, in a spin Hamiltonian through which 
equilibrium and out-of-equilibrium thermodynamics can 
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- in principle - be modeled. 

From a more applicative point of view, CoPhOMe 
could serve as a playground to better understand how 
thermal fluctuations affect broad DWs, typically hosted 
in ferromagnetic nanowires and described through a 
Hamiltoniari® 51 equivalent to Eq. ©• Domain walls 
in nanowires have been proposed as magnetic-memory 
elements to be man ipul ated by electric currentP^HSIl or 
injected spin waves 43 56 . Moreover, a detailed charac¬ 
terization of the electronic properties of these magnetic 
molecular helices is also relevant to further understand 
the interplay between structural chirality and magnetism 
in these materials. For instance , this interplay may give 
rise to inverse magneto-chiral effect 57 } i.e. magnetization 
induced by non-polarized light, or chiral-induced spin se¬ 
lectivity effects^ 59 , both of relevance for spintronics ap¬ 
plications. 
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Appendix A: Continuum limit 

Let us start considering the continuum limit for the M- 
NIT exchange interaction. Hamiltonian dla) is treated at 
the same level of approximation as Eq. T 2 ), namely Co 
spin operators are replaced by classical vectors |S 2 r | = 1 
and NIT spins by Pauli operators &2r+i m - 

d^M-NIT = — ^ [Jot $2 r ' & 2r+l + J (3 $2r ' **2r-l\ • 
r 

(Al) 

For a given configuration of the two classical spins (rep¬ 
resenting metal ions) S 2 r and S 2 r + 2 , the interlaying NIT 
spin experiences an exchange field , 1 = J a & 2 r + 
J( 3 ^ 2 r-\- 2 • The corresponding eigenvalues ±|h|^ +1 | de¬ 
pend parametrically on the orientation of both S 2 r and 
S 2 r+ 2 - Let us first compute the squared modulus of the 
local exchange field: 

|h^ +1 | 2 = (JcS 2r + J^2r + 2? 

= Jot + Jfi + ^JaJ/3^2r ' $2r+2 

= Jot + Jp + 2JaJ(3 ~ JotJ(3 (S 2 r — $ 2 r+ 2) 2 

— Ja + Jfi + 2J a Jp — J a J(3 {9 Z S) 2 

(A2) 


FIG. 9. Color online. x\\ — dM\\/dB\\ along the chain axis 
computed in B = 0 for 6 — 40°, 52°. Line-plus-symbol curves 
are obtained using the DFT values (set a) for J a ^ d/ 3 ; sym¬ 
bols correspond to calculations performed with a single, ef¬ 
fective value of J ex = 2JaJ/3/(Ja + d/ 3 ) (see the text). 


where we used the fact that classical spins have unitary 
modulus and substituted the difference between succes¬ 
sive classical spins by the spatial derivative. This approx¬ 
imation is essential to the continuum limit and relies on 
the assumption that large misalignment between S 2 r and 
S 2 r +2 is highly unlikely, because energetically not favor¬ 
able. Consistently, a Taylor expansion with respect to 
this derivative can also be performed to obtain 


ih^ +1 | = \j a + jf,wi- P (d z sy 


— | Ja T J3 | 


i-\ P (d z s ) 2 


(A3) 


with p = J a Jp/{Ja + d/ 3 ) 2 , always smaller than one. 
Small dev iations from the ground state of Hamilto¬ 
nian (Al) are described by the following functional of 
classical spins 


"Hex = / — p J ( d z s) 2 dz + const. (A4) 

Note that for J a = J(3 = d ex it is p = 1/4 and the 
exchange contribution given in Hamiltonian (|3| is ob¬ 
tained. Moreover, as long as the continuum approxima¬ 
tion is legitimate, one expects to obtain the same results 
using either J a ^ Jg or an effective unique coupling 
d e ff = 2 J a Jp/{Ja + d/ 3 ). We have seen that the set a 
of DFT calculations yields J a = —404 K and d /3 = —622 
K and that both are significantly larger than the esti¬ 
mate obtained for D z . Therefore, the continuum limit 
should be justified for the set of parameters considered 
in this work— 1 . If - in turn - all the excitations that are 
relevant for thermodynamics are contained in Hamilto¬ 
nian ©, it should be equivalent to compute, e.g., the 
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susceptibility using two different coupling J a = —404 K 
and Jq = —622 K or a unique effective coupling. This 
expectation is indeed confirmed quantitatively by Fig. [9] 
where curves computed using the two different values for 
J a and Jq given above or J e ^ = 2 J a Jp/( J a + Jp) — —490 
K are shown. This justifies why a single value of the ex¬ 
change coupling, rounded to J ex = —500 K, was used to 
study thermodynamic properties. 

In the following we show that in the continuum limit 
the anisotropy energy in Hamiltonian ([2]) transforms into 
the uniaxial anisotropy term of the functional |3|. Let 
us start rewriting the anisotropy energy as 

•Hani = -D z K • S ) 2 6(z-p- r)dz 

^ p,r 

= — D z f ^ {sin 2 9 [cos 2 <p(S x ) 2 + sin 2 p(S y ) 21 \ 

^ p,r 

-fcos 2 9(S Z ) 2 } S(z — p — r)dz 

(A5) 


where e r = (sin# cos (cp(r)) , sin#sin ((p(r)) , cos#) is the 
director of each local anisotropy axis (remember that 
c p = 27rr/3 imposed by the crystal symmetry) and 
S(. .. ) is the Dirac delta function, assuming unitary lat¬ 
tice units. As | J ex | is much larger than D z , it is rea¬ 
sonable to assume that the spin profile varies smoothly 
in low-energy configurations. Then, similarly to what 
is usual ly d one for spin spiral^, the cos 2 p and sin 2 p 
in Eq. (A5) can be replaced by their spatial averages, 
namely (cos 2 p} = (sin 2 ip)=1/2 (this is equivalent to as¬ 
suming a priori that DW excitations extend over several 
lattice units). The d epen dence on the site index (p + r) 
thus drops from Eq. (A5) as well as the summation and 
the Dirac delta function. More explicitly: 


n a ni = -D,f{\ Sin 2 0 [(S *) 2 + (S* 0 2 ] 
+ cos 2 9(S Z ) 2 } dz 
= -D z 

+ cos 2 #(S' Z ) 2 } dz 
= -D z 


i S in 2 #[l-(^) 2 ] 


J ^ sin 2 0 dz 


-D. 


/ 


cos 2 0 — ^ sin 2 0 


(S z ) 2 dz 


(A 6 ) 


where we have used the fact that ( S x ) 2 + ( S y ) 2 = 
1 — ( S z ) 2 . The second last line in Eq. ( |A 6 [ ) gives a 
contribution that does not depend on the orientation of 
the classical spins S(z); however, this term will gener¬ 
ally depend on the crystal structure, through 6. The 
last line gives the anisotropy contribution reported in 
Hamiltonian pi. The fact that the effective anisotropy 
D z (3 cos 2 9 — lj /2 is smaller than D z for any 9 ^ 0 and 
the agreement between the analytic model and transfer- 
matrix results reported in Fig. [ 8 ] justify a posteriori the 


validity of the approximations made to obtain Hamil¬ 
tonian © and the expression for the DW energy of 
CoPhOMe: £ dw (0) = ^2)^1 W z (3cos 2 6»-l). 


Appendix B: Transfer Matrix method 

The thermodynamic properties of classical-spin chains 
with nearest-neighbor interactions can be computed very 
efficie ntly with the (numerical) transfer-matrix tech¬ 
nique^^. This method can be extended to models in 
which classical and quantum spins alternate, like in our 
casd^. Noting that the quantum-spin operators are not 
directly coupled with each other, one can integrate out 
their degrees of freedom independently. Referring to 
Hamiltonian ([ 2 ]), we add the Zeeman contribution due 
to B to the exchange field acting on a generic quantum 
spin (introduced in Appendix 0 - which yields a total 
field hp 5 2 r+i = hp* 2 r+i + This field depends para¬ 

metrically on the orientation of the two classical spins 
S p? 2 r an d S Pj 2 r+ 2 - Note that the cell index p has been 
reintroduced, in accordance with Hamiltonian ©. The 
trace over the degrees of freedom of the quantum spin 
0 >, 2 r+i just brings a contribution 2 cosh (/3|h Pj 2 r +i|) into 
the partition function (here it is /3 = 1/T not to be con¬ 
fused with the label for the exchange coupling used in 
the rest of the paper). The total partition function reads 

z = J n p , r JHp 5 2r^C(Si 5 2, Si } 4) /C(Si } 4, Si } 6) 

X /C(S 1>6 , S 2 , 2 ) JC( S 2>2 , s 2>4 ) ... /C(S* P>6 , S 1>2 ), 

(Bl) 


where elementary kernels® are defined as 
^d(^p,2r? ^p,2r+2) — 2 COSh (/^|hp,2r+l I) 

X exp (^Sp i 2r-D2rSp,2r) X exp • Sp ? 2r) • 

(B2) 


Generally, the kernels /C(S Pj 2 , S Pj 4 ), /C(S Pj 4 , S Pj 6 ), 
/C(Si >6 , S 2 , 2 ) ar e not equivalent, which mirrors the fact 
that the three Co atoms in each cell are not equivalent 
from the magnetic point of view. This is crucial to model 
reciprocal non-collinearity of anisotropy axes. A kernel 
that involves only equivalent spins can, however, be ob¬ 
tained by tracing over the degrees of freedom of two Co 
spins in a generic cell, for instance those labeled with 4 
and 6 : 


^Tm(S Pj 2, S p+ 1^2) — J dfl P: 4 J /C(Sp 5 2 , S p> 4) 

x /C(S p? 4, S p,o) A 5 (S p?6 , S p _|_i,2) dQpfi . 


(B3) 


Then, the calculation of the partition function in 
Eq. (Bl) reduces to the following eigenvalue problem 


/ 


JC TM (Si, S i+1 )W m (S i+1 )dtt i+1 = X m W m (Si ), (B4) 
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whose eigenvalues may typically be ordered from the 
largest to the smallest one: Ao > Aj > A 2 > .... The 
integrals on the solid angles were performed numerically 
by discretizing the u nitar y sphere^HMl As a result, the 
eigenvalue problem ( |B4| ) was converted into a linear- 
algebra problem 67 . The resulting matrix - to be diag¬ 
onalized - is generally not symmetric and its eigenvalues 
are pairs of complex conjugates. In the thermodynamic 
limit N p 00 the partition function equals the N p - th 
power of the largest eigenvalue, Z = A 0 p (N p is the 
number of cells); therefore Ao has to be real. For each 
temperature, the number of points used to sample solid 
angles was increased until the free-energy F = —Tin (Z) 
converged to a stable value 4 -R From the free-energy, the 
magnetization per unit cell M p = —dF/ ( N p dB p ) and the 
susceptibility Xv — dM^/dB^ were obtained (77 = x,y,z 
indicates the crystallographic axis along which B is ap¬ 
plied). 


order of the correlation length within a lapse of time 
equal to the relaxation time, one gets to the equiv¬ 
alence £ 2 ~ Dr. As already reported in Ref. 06J 
for the Ising limit corresponding to 5 = 1 , this map¬ 
ping yields the same asymptotic behavior obtained by 
Glauber: = 2^ 2 /a. The subscript indicates that this 

result actually holds for the infinite chain. For a finite 
chain, still in the framework of Glauber dynamics, it can 
be shown that the low-temperature behavior of the re¬ 
laxation time is r ~ Le Sd ™' T /(2a), L being the chain 
length 6 ^. Our goal is now to estimate what effective a 
should be used to mimic the diffusion of a broad DW with 
Glauber dynamics. The first step is to identify S with the 
DW width w\ for the Hamiltonian ©> , the latter reads 
w = VP37 [21^(3cos 2 6 — 1)]. Adap ting the numerical 
factors given in Ref. Hi and 02] to Eq. ( |C3| ), the diffusion 
coefficient can be related to the spin-Hamiltonian param¬ 
eters: 


Appendix C: Glauber model and diffusion of broad 
domain walls 


D ~ 0.34 


w 2 T 
Td £ dw ( 0 ) 


(C4) 


One of the reasons for the fortune of the Glauber model 
in the SCM community certainly lies on the possibility of 
deriving analytic results. On the contrary, the treatment 
of stochastic dynamics in the framework of the Heisen¬ 
berg model is not as simple and even numerical simula¬ 
tions - involving a physically meaningful time scale - are 
challenging. The starting point is the stochastic Landau- 
Lifshitz-Gilbert equation, which can be integrated fol¬ 
lowing the Langevin or the equivalent time-quantified 
Monte-Carlo approach ! 42 1 68 1 . Here we limit ourselves to 
propose a mapping of the elementary time scale of the 
Glauber model into the diffusion coefficient associated 
with the motion of a broad DW. Within the Glauber 
model the probability to flip a spin occupying the site i 
of an Ising chain is given by 21 


a 

Wi = — 
2 


1 — 2 a i { a i +1 + a i-1) 


(Cl) 


where cq = ±1 are Ising spin variables and 7 = 
tanh (2 J ex /T). When the i-th spin is the front edge of a 
DW, it is cq + 1 + cq_i = 0; therefore, a/2 can be identi¬ 
fied with the probability to move a DW of one site per 
unit time. Let us define P(z,t) as the probability to 
find a DW at the coordinate 2 - being now a continuous 
variable - at time t. This probability obeys the master 
equation 

d t p(z, t) » I [P(z + S,t) + P{z - 5, t) - 2 P(z, t)] (C2) 


where 5 is the elementary d ispl acement of the stochastic 
process. By expanding Eq. (C2) for small S, the standard 
ID diffusion equation is recovered: 


d t P(z,t) = ^S 2 d 2 P(z,t) (C3) 

in which D = ad 2 /2 plays the role of diffusion coeffi¬ 
cient. Assuming that a DW covers a distance of the 


where Td is the damping time ra ~ Ti/{2D z ao), ao is 
the Gilbert damping 70 'and D z is the anisotropy on each 
CojUH] For this special case the attempt frequency of 
the Glauber model and the diffusion coefficient are re¬ 
lated by the equation a = 2 D/w 2 . The low-temperature 
expansion given above for the relaxation time of a finite 
Ising chain can tentatively be extended to a chain with 
broad DWs: 

w 2a 2 w 0.68 2 V™ 0.68/ 5 V J 


in the last passage the approximation ^e K ^ 6 e 1,06/ ^ was 
made based on the fact that the experimentally relevant 
range of k = £dw(0)/T is 12 < n < 22. In Eq. (C5), 
the inverse of the term in parenthesis represents the ef¬ 
fective Glauber attempt frequency that should emulate 
the diffusion of a broad DW. For our computation pa¬ 
rameters and 50° < 0 < 54°, it is a e ff ~ 10 15 <aG s _1 . For 
the temperature range where slow relaxation is observed 
in CoPhOMe, the order of magnitude of the Gilbert 
damping can be roughly estimated from EPR spectra 
of the C 0 NIT 2 , shown in Ref. [39| In a trivial preces- 
sional model, the half-with-half-maximum AH of the res¬ 
onance peak would be related to the Gilbert damping as 
7 g AiL ~ with c^l = 7 g^o being the Larmor fre¬ 

quency, Hq the resonance field and y g the gyromagnetic 
ratio. Then one obtains gq — AH/Ho — 0.15, from 
which a e ff ^ 1.5-10 14 s -1 . To simulate the data reported 
in Ref. 17 the Glauber model with a = 2.6 • 10 13 s -1 was 
used. Considering the crude approximation made to es¬ 
timate gg and that the parameters in Hamiltonian ([2| 
have been deduced from DFT calculations and not fit¬ 
ted, the present mapping sounds in reasonable agreement 
with previously published results. 
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